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This article presents a system dynamics model to analyze the growth of cars 
and the effect of different policies on carbon emissions from the transport 
sector. The simulation model used in this work was built using the 
methodology of systems dynamics (SD) developed by Jay W. Forrester at the 
Massachusetts Institute of Technology (MIT). The model was applied to the 
transport sector of the city of Bogota, Colombia for a period of time between 
2005 and 2050. The information used to feed the model comes from reliable 
sources such as DANE (National Administrative Department of Statistics) 
and EIA (U.S Energy Information Administration). Four scenarios were 
proposed that relate urban development policy and environmental policy. 
The main results indicate that the number of cars in Bogota can reach up to 


13 million vehicles in 2050 and the projection of CO2 emissions would reach 


Transport sector model tee h 1 i 
34 million TonCO2 in the absence of an appropriate environmental policy. 


Transport simulation 
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1. INTRODUCTION 

Transport plays an important role in the social and economic development of a country. But on the 
other hand, in relation to the environment it is a source of emissions, noise and vibration and causes risks to 
health and safety [1]. As the sector with the highest oil consumption and the fastest oil demand and carbon 
emissions growth, transport needs to improve energy efficiency and the development of low carbon 
technologies [2]. 

Transport allows great advantages to people in terms of commuting, but also consumes a large 
amount of energy and emits great amounts of greenhouse gases (GHGs) and pollutants. Because cities are the 
centers of human transport, urban transport has become the sector with the highest energy consumption and 
GHG emissions generation [3]. The United Nations reported that a fifth of the world's population lives in 600 
cities in the world. It is expected that the population of cities in developing countries will double by 2030 and 
that their energy consumption will represent 60% to 80% of the total world energy consumption. Thus, GHG 
emissions will represent 70% of the total human GHG emissions with main emission sources concentrated in 
fossil fuels for energy supply and the transport sector [4]. 

Many academics have carried out studies to predict energy demand and carbon emissions. 
According to the different research objectives, these can be divided into three categories [4]: First, there are 
predictions of energy consumption and carbon emissions at the global, regional and national levels, as is the 
case with the annual global energy and carbon emissions of the International Energy Agency (IEA) [5]. 
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Second, there are predictions of energy consumption and carbon emissions in specific sectors of different 
regions in certain countries [6]-[8]. In this sense, [9] combines the time-heat-power optimization model and 
the optimization model of a decentralized energy system to forecast energy consumption in the construction 
sector in Germany. In [10] the historical trend of the energy demand of road transport in developed 
economies is analyzed and the trajectory analysis method and the BMA model (Bayesian Model Averaging) 
were applied to analyze the energy demand of road transport. In [11] Autoregressive—moving-average model 
(ARMA) was used to forecast the energy demand and environmental emissions of India's iron manufacturing 
industry. In [12] a linear logarithmic equation was applied to explain the relationship between carbon dioxide 
emissions and their influence factors, to forecast the CO2 emissions of China's electric power industry during 
the 2016-2030 periods under different scenarios. Third, there are predictions of specific kinds of energy at 
different scales. For example, in [13] China's demand for natural gas is predicted during the period 2015- 
2020 using a self-adapting intelligent grey prediction model. In [14] the future evolution of the energy 
structure of Brazil is analyzed and the demand and supply of power is predicted to 2030. 

There are researches that apply different methods and models [15]-[25] to evaluate the demand and 
energy consumption of the transport sector of different cities and countries. The methodologies allow 
predicting future transport emissions as well as energy demand under different possible scenarios. This paper 
presents a model of systems dynamics using different scenarios to analyze two variables of the transport 
sector of the city of Bogota: automotive fleet growth and emissions generation. 


2. MODEL DESCRIPTION 

The simulation model used in this work was built using the methodology of systems dynamics (SD), 
this methodology was developed by Jay W. Forrester at the Massachusetts Institute of Technology (MIT) 
around the year 1950, SD is a tool that employs a system of differential equations that allows to establish and 
understand the behavior of complex systems in order to propose future scenarios and make decisions [26]. 

SD is appropriate for modeling social systems under conditions of uncertainty and imbalance [27]. 
Cities are social systems that present such conditions, since they face uncertainty associated with urban 
growth and regulatory interventions at the environmental level [27]. Consequently, the SD is an adequate tool 
to analyze the effect of this uncertainty on the cycles of growth and planning in the long term. In particular, 
the simulation model proposed constitutes a laboratory for policy analysis, since it studies how urban growth 
and the transport sector react to policy changes [26], [28]. 

The problem of emissions generation of the transport sector can be seen as a diagram of subsystems. 
The diagrams of subsystems facilitate the definition of the limits and aggregation degree of the model, by 
allowing to define the structure of the model as a whole composed by modules, where each module 
represents a physical or conceptual unit relevant for decision making [26]. Next, the modules shown in 
Figure 1 are described. 
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Figure 1. Structure of the simulation model represented in a subsystem diagram 
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The population module depends on the entry (immigration, births) and departure (deaths, migration) 
of the population. The vehicles module depends on the number of passengers and the number of trips they 
make. The emissions module depends on the total number of trips and the type of vehicle in which these trips 
are made, since each type of vehicle has a different emission factor. Finally, the passenger’s module depends 
on the amount of population that is mobilized and the number of trips made during the year. 

The stocks and flows diagram of the model is presented in Figure 2, which also shows an overview 
of the feedback loops. Following, the investment cycles B1, B2, R1 and R2 will be described. 
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Figure 2. Stocks and flows diagram 


The cycle of private vehicles B1 shows the behavior of the purchase of these vehicles. An increase 
in the number of private vehicles will produce an increase in traffic density, which in turn makes the 
purchase of private vehicles less attractive. This last effect will produce a reduction in the purchase of private 
vehicles. 

Cycles R1, R2 and B2 correspond to the purchase of vehicles for public transport. In cycle R1, an 
increase in the number of public transport vehicles generates a greater sense of comfort among users (this is 
represented by the variable "demand vehicle effect"), which in turn increases the journeys of people and 
therefore the number of vehicles required. This last effect increases the purchase of public transport vehicles. 
In cycle R2, an increase in the number of public transport vehicles generates a higher density of traffic, which 
increases the trips of people and therefore the number of vehicles required, the latter effect increases the 
purchase of public transport vehicles. In cycle B2, an increase in the number of public transport vehicles 
reduces the number of vehicles required, which reduces the purchase of public transport vehicles. The main 
equations of the model are described in Table 1. 
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Table 1. Equations of the Simulation Model 


Equation Units 
Line felt) + DT,(t) (1) 


The emissions of the transport sector depand on the distance traveled by type of technology Million tons CO2 
(DT) and the emission factor of each technology (f). E T: 


DT;(t) = DTT,(t)+ NT,(t)+NTV,(t) (2) 
The total distance traveled (DT) by type of vehicle depends on the number of trips made by Ka 
each vehicle on a trip (DTT), the total number of trips made and the total number of vehicles. 


NTV,(t) = J NVR,(t) dt G) 
The total number of vehicles (NTV) depands on the purchase of vehicles or number of vehicles 


iced (NVR). Vehicles 
NVR, (t) = MAX (0, NP,(T) + UtilMean(t) — NVR; (t— 1)) (4) 
The number of vehicles required at time t (NVR) depends on the number of trips (NP), and the Vehicles 
capacity of the vehicles according to the type. 
Po(t) = {Change(t)+Po(t—i)dt (5) 
The population (Po) depands on the percentage of change (Change) and the population in the Person 


previous period. 


The information used to feed the model comes from reliable sources such as DANE (National 
Administrative Department of Statistics) and EIA (U.S Energy Information Administration). The simulation 
period comprises from 2005 to 2050. In the model, private and public vehicles are considered. Private 
vehicles are divided into: cars, motorcycles and electric bicycles. Public vehicles are divided into: BRT (Bus 
Rapid Transit, for the study case of Transmilenio), buses, taxis and metro. In the model it is considered that 
the metro will begin operations in year 2020 and then every 4 years a line will enter according to the scenario 
(See section 3). In Table 2, the main parameters used and their corresponding values are shown. 


Table 2. Values of Parameters 


Parameters Value Units Ref. 
Initial Population 6.840.120 Person 29] 
Initial private vehicles 698.893 Vehicles 30] 
Avg gasoline vehicle fuel efficiency 11/100 L/Km 31] 
Avg diesel vehicle fuel efficiency 9,8/100 L/Km 31] 
Avg natural gas vehicle fuel efficiency 0.011/100 m3/km 31] 
Avg diesel BRT fuel efficiency 35/100 L/Km 31] 
Avg diesel bus fuel efficiency 25/100 L/Km 31] 
Avg gasoline Taxi fuel efficiency 9/100 L/Km 31] 
Avg natural gas Taxi fuel efficiency 0.011/100 m3/Km 31] 
Energy content of gasoline 0.0344 GJ/L 32] 
Energy content of diesel 0.0371 GJ/L 32] 
Energy content of Natural gas 0.0356 GJ/m? 32] 
Gasoline emission factor 69,25/1000  tonCO2/GJ 33] 
Diesel emission factor 74,01/1000  tonC0O2/GJ 33] 
Natural gas emission factor 56,1/1000 tonCO2/GJ 33] 


The following section presents the results of the simulation model that correspond to the emissions 
of the transport sector in the city of Bogota under different scenarios. 


3. RESULTS 

Four simulation scenarios have been proposed, which are the result of several workshops conducted 
in Colombia with mobility experts. The four scenarios presented in Figure 3 combine the two most uncertain 
variables in the system: urban development policy and environmental policy, which are represented on the x 
and y axis. The urban development policy encourages construction on the territory. Environmental policy 
promotes cleaner forms of transport. 

Scenario 1 is the optimistic scenario, in which an urban and environmental development policy is 
applied. Under this scenario, 55% of urbanization and 45% of public and green areas are considered, in 
addition to 5 metro lines. 
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Scenario 2 gives greater importance to urban development policy and does not consider an 
environmental policy. In this scenario, 80% of urbanization and 20% of public and green areas are 
considered; additionally the metro will have 2 lines. 

Scenario 3 has an environmental policy. In this scenario only 30% is urbanized and 70% is delimited 
for public and green areas, in addition to 3 metro lines. 

Scenario 4 is the pessimistic scenario. Under this scenario high urbanization of the city is considered 
corresponding to 90% of the available territory and little availability of green areas, in addition to 1 metro 
line. 

On these scenarios the number of private vehicles, the total emissions of the transport sector and the 
contribution of this sector on the total emissions of the city have been analyzed. 
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Figure 3. Scenarios to study the emissions of the transport sector under different levels of urban growth 


14 35 
12 30 
= 
a 10 
Q 25 
2 8 
£ 
g Q 20 
6 [0] 
Oo 
4 5 15 
FE 
2 s 10 
0 5 
` N 
SSSSSSSSSAASSSRRSSRTTITTE EEE ECE EEEEEELELELTELE 
REE Ae a Ren Se Ae gy ce ewe OP Pee RRIRRRRRRRRRRRRRRRRRRRRRR 
—— scenariol = scenario? = scenario3 == scenario 4 = scenariol = scenario2 === scenario3 === scenario 4 


Figure 4. Growth of the automotive fleet in Bogota Figure 5. Total emissions from the transport sector in 


Bogota 

85% 

80% 

75% 

70% 

65% 

60% 
woOnRaARATOWON TO WA MNRA A MNHWAOS 
oo Odd at ANNAN A AM MMMM FFAN 
soaooo0oo0oooo0oooooooooo coco ood 
NANN NN NN NN N N NN N N N NN NANNAN 

—— scenariol == scenario2 === scenario3 =< scenario 4 


Figure 6. Participation of the transport sector in the total emissions of the city of Bogota 
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Figure 4 shows the automotive fleet growth in the city of Bogota. It is possible to observe that in 
scenarios 2, 3 and 4 there is a great growth in the number of vehicles due to the lack of an environmental 
policy. 

On the other hand, in scenario 1 ("optimistic" scenario) the presence of an environmental policy 
leads to greater use of public transport and therefore fewer vehicles. Figure 5 shows the growth of carbon 
emissions in the transport sector in Bogota. Because scenarios 2 and 4 do not have an environmental policy, 
there is a faster growth of emissions due to the tise on the vehicle fleet and the lower participation of the 
metro. 

In contrast, scenarios 1 and 3 do consider an environmental policy. Under these scenarios, in 
scenario | emissions grow less than in scenario 3. In addition, scenario | shows growth at the beginning of 
the simulation period and then a decrease in emissions. This decline is due to a greater participation of the 
metro which is a clean mean of transport. Figure 6 shows the participation of the transport sector in the total 
emissions of the city of Bogota. In scenarios 2 and 4 an environmental policy is not considered, therefore the 
transport sector is responsible for most of the total emissions of the city of Bogota. Under scenario 3, 
emissions from the transport sector show a growing trend. Finally, in scenario 1 (optimistic scenario), most 
of the emissions come from the transport sector, although the percentage of participation is lower than in the 
other scenarios. 


4. CONCLUSIONS 

This article uses an SD model to learn about the design of policies aimed to reduce emissions in the 
transport sector. Several lessons are derived from this research from the point of view of sustainable planning 
of cities. In particular, a comprehensive urbanization policy is necessary to effectively reduce emissions from 
the transport sector. This comprehensive policy implies the coexistence of a policy for urban development 
and an environmental policy. 

The urban development policy facilitates the controlled growth of the city, keeping emissions from 
the transport sector under acceptable levels. Additionally, this policy is a complement to the environmental 
policy. The key variable in the design of an environmental policy is the construction of a metro in the city of 
Bogota, because this would significantly reduce carbon emissions, as indicated by the simulation results. 

The SD model built for this research complies with the objective for which it was designed, that is, it 
allows to analyze the growth of the vehicle fleet and the effect of different policies on the carbon emissions 
of the transport sector. The control of carbon emissions from the transport sector is crucial, since this is the 
sector with the highest share of total emissions. Therefore, it is necessary to study in the future the design of 
clear policies that help to encourage, other types of mobility, such as electric mobility and the use of bicycles. 
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